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Conditions  of  Use 


10  December  1984 


The  following  conditions  regulate  the  use  of  computer  programs 
developed  by  the  Hydrologic  Engineering  Center  (HEC),  Corps  of 
Engineers,  Department  of  the  Army. 

1 .  The  computer  programs  are  furnished  by  the  Government  and  are 
accepted  and  used  by  the  recipient  individual  or  group  entity 
with  the  express  understanding  that  the  United  States  Government 
makes  no  warranties,  expressed  or  implied,  concerning  the 
accuracy,  completeness,  reliability,  usability,  or  suitability 
for  any  particular  purpose  of  the  information  or  data  contained 
in  the  programs,  or  furnished  in  connection  therewith,  and  that 
the  United  States  Government  shall  be  under  no  liability 
whatsoever  to  any  individual  or  group  entity  by  reason  of  any  use 
made  thereof . 

2.  The  programs  belong  to  the  United  States  Government. 
Therefore,  the  recipient  agrees  not  to  assert  any  proprietary 
rights  thereto  nor  to  represent  the  programs  to  anyone  as  other 
than  Government  programs . 

3.  The  recipient  may  impose  fees  on  clients  only  for  ordinary 
charges  for  applying  and  modifying  these  programs. 

4.  should  the  recipient  make  any  modifications  to  the 
program(s),  the  HEC  must  be  informed  as  to  the  nature  and  extent 
of  those  modifications.  Recipients  who  modify  HEC  computer 
programs  assume  all  responsibility  for  problems  arising  from,  or 
related  to,  those  modifications.  User  support  from  the  HEC  to 
third  party  recipients  will  only  be  provided  after  the  recond 
party  demonstrates  that  program  difficulties  were  not  caused  by 
their  modifications. 

5.  This  "Conditions  of  Use"  statement  shall  be  furnished  to  all 
third  parties  that  receive  copies  of  HEC  programs  from  the 
recipient.  Third  party  recipients  must  be  notified  that  they 
will  not  receive  routine  program  updates,  correction  notices,  and 
other  program  services  from  the  HEC  unless  they  obtain  the 
program(s)  directly  from  the  HEC. 

6.  All  documents  and  reports  conveying  information  obtained  as  a 
result  of  the  use  of  the  program(s)  by  the  recipient,  or  others, 
will  acknowledge  the  Hydrologic  Engineering  Center,  Corps  of 
Engineers,  Department  of  the  Army,  as  the  origin  of  the 
program (a ) . 
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COMPUTER  IMPLEMENTATION  GUIDE  FOR  HEC-1 


1 . 0  INTRODUCTION 

1.1  Purpose  of  Document 


This  guide  provides  information  about  the  implementation,  organization  and 
structure  of  the  source  code  for  the  HEC-1  Flood  Hydrograph  Package,  1981.* 

It  is  intended  for  use  by  systems  engineers  and  engineering  programmers  who 
are  implementing  the  program  on  a  mainframe  computer  system  or  modifying  the 
source  code  to  meet  special  needs  of  their  users  and/or  computer  system. 

Some  of  the  information  in  this  Implementation  Guide  will  be  useful  to 
persons  modifying  the  microcomputer  (PC)  version  of  HEC-1.  Separate  documents 
are  available  for  installation  of  the  PC  executable  code,  "Microcomputer 
Version  of  HEC-1  Flood  Hydrograph  Package,"  and  for  use  of  the  PC  source  code, 
"Source  Code  Comments:  HEC-1  Microcomputer  Version.”' 

1.2  Supplementary  ProRrams 

A  program  to  convert  old  input  data  (1973  version)  into  the  current  input 
data  format  (1981  version)  is  available.  Corps  users  may  interface  HEC-1 
input  and  output  with  the  HEC  Data  Storage  System.  These  programs  are 
described  in  the  Program  Implementation  section  that  follows. 

1.3  Notification  of  Errors  and  Modifications 


All  holders  and  users  of  the  HEC-1  program  will  be  notified  of  errors  in 
or  modifications  to  the  program.  The  HEC  maintains  an  address  file  of  persons 
who  have  obtained  the  program  from  HEC  as  well  as  those  people  who  use  the 
program  on  a  computer  service  and  have  made  official  arrangements  with  the  HEC 
to  be  placed  on  the  HEC-1  mailing  list. 

Users  of  the  program  should  notify  the  HEC  about  any  errors  they  find  in 
the  program.  The  HEC  also  welcomes  any  suggestions  for  improvement  of  the 
program  source  code,  operating  mode  and  technical  capabilities.  Any  informa¬ 
tion  regarding  problems  with  the  use  of  the  program  on  different  computers 
will  be  greatly  appreciated. 

2.0  PROGRAM  IMPLEMENTATION 


2.1  Using  the  HEC-Supplied  Magnetic  Tape 


The  HEC  distributes  program  source  code  for  mainframe  computers  on  a  1/2" 
magnetic  tape  written  according  to  the  requester's  specifications  (for  a  PC,  a 
5-1/4"  floppy  disk).  A  description  of  the  tape  contents  Is  provided  with  the 
tape.  The  files  on  the  tape  are:  source  code,  test  data,  input  data 
description,  and  a  data  conversion  program. 


*or 


IP' 


•>a 


*HEC-1  Flood  Hydrograph  Package,  Users  Manual,  U.S.  Army  Corps  of  Engineers,  o/ 


Davis,  California,  September  1981,  Revised  January  1985. 


**ty  Cedes 


Diet 


Avail  and/or 
Spool*! 


OO 


The  source  code  should  be  copied  from  the  tape  into  the  user's  program 
library.  The  source  code  has  standard  FORTRAN  77  headings  for  the  main  routine 
and  subroutines.  Some  editing  of  the  source  code  may  be  necessary  as  described 
in  the  next  section  on  machine-dependent  Code  and  Memory  Storage  Requirements. 
The  Input  Data  Conversion  Program  is  a  separate  program  and  should  be  compiled 
separately  from  the  HEC-1  source  coda. 

2.2  Input/Output  File  Structure 

HEC-1  uses  up  to  14  input/output  and  scratch  files.  These  can  be  stored 
on  disk,  tape,  or  whatever  medium  is  available.  The  unit  numbers  assigned  to 
the  HEC-1  files  are  shown  in  Table  1.  These  numbers  can  be  changed  for  a 
particular  installation  by  changing  their  definition  in  BLOCK  DATA.  Table  II 
describes  the  execution  time  file  assignments  for  several  systems. 

2 . 3  Machine-Dependent  Code 

Two  items  in  the  source  code  are  unique  to  various  computer  systems  as 
shown  in  Table  III.  The  first  item  is  necessary  when  using  FORTRAN  77.  The 
OPEN  statement  shown  specifies  that  a  blank  in  the  HEC-1  input  data  is  read  as 
a  zero.  The  second  item  is  a  system-related  feature  used  in  subroutine  BANNER 
to  obtain  the  date  and  time  for  a  particular  program  execution.  The  subrou¬ 
tines  DATE  and  TIME  shown  are  examples  of  this  feature  on  a  CDC  Cyber  865 
computer  system. 

2.4  Compilation  and  Execution  Requirements 

HEC-1  requires  a  Fortran  77  compiler  and  up  to  14  input/output  files.  The 
computer  memory  requirements,  compilation  time  and  example  problems'  execution 
times  are  given  in  the  HEC-1  Users  Manual  (Section  13,  Table  13.2).  The  HEC 
would  appreciate  information  about  memory  requirements  and  execution  times  on 
other  computers. 

2.5  Program  Operation  Verification 

The  correct  operation  of  HEC-1  should  be  verified  using  the  HEC-supplied 
test  data.  The  te?*-  data  corresponds  to  the  example  problems  provided  in 
Section  12  of  the  HEC-1  Users  Manual.  The  results  reported  for  the  first  11 
test  examples  should  be  reproduced  exactly,  with  minor  exceptions  due  to  word 
size.  Test  example  12  is  more  sensitive  to  word  size.  Comparison  with  the 
reported  results  should  be  within  five  percent.  The  reported  results  were 
produced  on  the  Harris  500  computer  using  a  39-bit  mantissa  for  real  numbers. 

2.6  Input  Data  Conversion  Program 

The  Input  Data  Conversion  program  provided  on  the  tape  converts  old  HEC-1 
(1973)  data  sets  into  the  equivalent  data  sets  required  by  the  new  (1981) 
program.  The  conversion  program  does  not  recognize  the  special  records  used 
in  the  old  dam-break,  kinematic-wave  and  system-optimization  versions  of  HEC-1. 
Input  to  the  conversion  program  is  the  old  data  set  (A,  B,  etc.,  records)  and 
the  new  data  set  is  output  as  a  scratch  file  for  disposal  to  the  card  punch, 
permanent  file  and/or  printer.  The  conversion  program  has  not  been  extensively 
tested  but  should  suffice  to  make  the  majority  of  the  changes  necessary  to  use 
the  new  program.  Users  are  encouraged  to  make  the  conversion  and  then  discard 
the  old  data  deck.  The  conversion  program  is  not  intended  to  be  a  preprocessor 
to  the  new  program. 


2 


Table  I 


Input/Output  and  Scratch  Files 


Unit 

Number 

Variable 

Name 

Description 

Formatted,  F 
Unformatted.  U 

Max  Record  Lenqth 

5 

INP 

Primary  input 

F 

80  characters 

6 

IP 

Primary  output  file  (printer) 

F 

132  characters 

7 

IPU 

Punch 

F 

80  characters 

23 

IC 

Working  input  file;  reformatted  input  data  with  line 
number  and  next  record  ID  appended  to  front  of  each 

record 

F 

89  characters 

24 

IS* ** 

Dam-overtopping  summary  report 

F 

132  characters 

25 

IU* 

Runoff  parameter  optimization 

F 

132  characters 

32 

IDIV 

Scratch,  saves  diversion  hydrographs 

U 

4895  real  -*■  3  integer  words 

33 

IE 

Scratch;  expected  annual  damage  summary  data 

U 

50  real  +  6  integer  words 

34 

IR 

Scratch;  data  for  first  plan  in  multiplan  run 

u 

61  real  words 

35 

I  SOP 

Scratch;  data  for  flood  control  system  optimization 

u 

2400  real  words 

36 

LSFIL 

Scratch;  data  for  user-defined  output  tables 

u 

301  real  words 

38 

ND 

Scratch;  output  summary  data 

u 

91  real  +  4  integer  words 

** 

IOUT 

Output  data;  used  to  save  hydrographs  for  a  subsequent  job 

F 

131  characters 

** 

IQIN 

Input  data;  hydrographs  from  a  previous  job 

F 

131  characters 

*  File  is  copied  to  primary  output  file  (IP)  by  subroutine  PRT 

**  Unit  number  is  defined  by  user  on  KO  or  BI  records  (The  unit  numbers  specified  should  not 
conflict  with  other  file  definitions,  for  example,  21  and  22  are  possible  choices). 
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Table  II 


Input/Output  File  Specifications  for 
Alternative  Computer  Systems 


Table  Ila.  File  Assignments  for  Using  HEC-1  on  an  IBM  370 


//HEC1N  PROC  PR0GRAM=HEC1N , REG=512K ,LIB=TESTLIB 

//HEC1N  EXEC  PGM=&PR0GRAM,REGI0N=&REG 

//STEPLIB  DD  DSN=WATER . &LIB , DISP=SHR 

//FT01F001  DD  DDNAME=SYSIN 

//FT02F001  DD  SYS0UT=B 

//FT03F001  DD  SYSOUT=A,DCB=(LRECL=133 , RECFM=FA, BLKSIZE=133 ) 

//FT10F001  DD  UNIT=DISK,DCB=(RECFM=FB , LRECL=89 , BLKSIZE=3916 ) , 
//  SPACE= ( CYL , ( 2 , 1 } ) 

//FT11F001  DD  UNIT=DISK , SPACE=( CYL , { 1 , 1 ) ) , DCB=RECFM=VBS 

/ /FT12F001  DD  UNIT=DISK , SPACE= ( CYL ,(1,1)), DCB=RECFM=VBS 

//FT14F001  DD  UNIT=DIS.K , SPACE=  ( CYL ,  ( 1 , 1 ) )  , DCB=RECFM=VBS 

/ /FT16F001  DD  UNIT=DISK,SPACE=(TRK, (10, 1) )  , 

DCB=(RECFM=FB , LRECL=132 , BLKSIZE=3960) 

//FT17F001  DD  UNIT=DISK ,SPACE=(CYL, (1,1)), DCB=RECFM=VBS 

//FT18F001  DD  UNIT=DISK,SPACE=(CYL, (1 ,1) > , DCB=RECFM=VBS 

//FT19F001  DD  UNIT=DISK , SPACE= (TRK, ( 10 , 1 ) ) , 

DCB= ( RECFM=FB , LRECL=132 , BLKSIZE=3960 ) 

//FT20F001  DD  UNIT=DISK , SPACE=(CYL , ( 1 , 1 ) ) , DCB=RECFM=VBS 


UNIT  (FILE)  NUMBERS  ARE  IMPLEMENTATION  DEPENDENT  AND  WILL  VARY  BETWEEN  SITES. 
FOR  THIS  EXAMPLE  THE  UNIT  NUMBERS  AS  ASSIGNED  IN  BLOCK  DATA  ARE: 

INP  IP  IPU  IC  IS  IU  IDIV  IE  IR  ISOP  LSFIL  ND 

1  3  2  10  16  19  11  18  17  20  14  12 


Table  II. b  File  Assignments  for  Using  HEC-1  on  a  Burroughs  B6700 

FILE  5  ( KIND=DISK ,  FILETYPE=7  , TITLE="D/HEC1” ) 

FILE  6=FILE6 , UNIT=PRINTER , REC0RD=22 

FILE  7=FILE7 ,UNIT=PRINTER,RECORD=22 

FILE  11=FILE11 ,  UNIT=DISK ,  RECORD=22  ,  BI.OCKING=30 

FILE  12 -FILE12 ,UNIT=DISK , REC0RD=22 , BLOCKING=30 , AREA=120*500 

FILE  23=FILE23 ,UNIT=DISK, REC0RD=15 , BLOCKING=30 , AREA=120*500 

FILE  24=FILE24 , UNIT=DISK , RECORD=22 , BLOCKING=30 , AREA=120*500 

FILE  25=FILE25 ,UNIT=DISK , REC0RD=22 , BLOCKING=30,AREA=120*500 

FILE  32=FILE32 ,UNIT=DISK , REC0RD=4895 , AREA=1*1000 

FILE  33=FILE33 ,UNIT=DISK , REC0RD=60 , BLOCKING=30 , AREA=120*500 

FILE  34=FILE34 ,UNIT=DISK ,RECORD=90 ,BLOCKING=30 ,AREA=1 20*500 

FILE  35=FILE35 ,UNIT=DISK , REC0RD=2400,AREA= 1*1000 

FILE  36=FILE36 ,UNIT=DISK, RECORD=330 ,BL0CKING=6 , AREA=120*500 

FILE  38=FILE38 ,UNIT=DISK,REC0RD=120 , BLOCKING=30 , AREA=120*500 
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Table  lie.  File  Assignment  for  Using  HEC-1  on  a  CPC 


Place  the  following  lines  as  the  first  record  of  FORTRAN  program  source  code: 

PROGRAM  HEC1  ( INPUT .OUTPUT , PUNCH , TAPE5=INP0T , 

1  TAPE6=OUTPUT,TAPE7=PUNCH,TAPE21 .TAPE22 .TAPE23 .TAPE2A , 

2  TAPE25 .TAPE32 .TAPE33 .TAPE34 .TAPE35 .TAPE36 , 

3  TAPE38 ) 


See  Table  I  for  unit  numbers  corresponding  to  TAPES  through  TAPE38.  For 
example,  TAPE5  corresponds  to  unit  number  5,  variable  INP,  and  TAPE38 
corresponds  to  unit  number  38,  variable  ND. 


Table  lid. 

File  Assignments 

for  Using  HEC-1  on  a  Harris  500 

AS 

5  =  Input  File 

formatted  file 

AS 

6  =  Output  File 

M 

AS 

7  =  Punch  File 

M 

AS 

23  =  W1 

formatted  work  file 

AS 

2A  =  W2 

t# 

AS 

25  =  W3 

•t 

AS 

32  =  U2 

unformatted  work  file 

AS 

33  =  U3 

M 

AS 

3A  =  UA 

•• 

AS 

35  =  U5 

ft 

AS 

36  =  U6 

•« 

AS 

38  =  U8 

M 
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Table  III 


Machine-Dependent  Code 


PROGRAM  HEC1 


C  USE  FOLLOWING  STATEMENT  IF  PROGRAM  IS  COMPILED  USING 

C  A  FORTRAN  77  COMPILER. 

C 

OPEN  (IC,BLANK=* ZERO') 

C 


1147002 

1147102 

1147202 

1147302 

1147402 

1147502 

1147602 


SUBROUTINE  BANNER 


1480 

MACHINE  SPECIFIC  CALL  TO  DATE  AND  TIME  ROUTINES  1481 

1482 

CALL  DATE  (II)  1483 

CALL  TIME  (12)  1484 
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2 . 7  Data  Storage  System  -  DSS 


The  HEC  Data  Storage  System  is  not  currently  available  for  general  distri¬ 
bution.  Consequently,  an  HEC-1  dummy  routine,  DSMSTR,  is  used  in  place  of  the 
routines  available  to  Corps  users  for  accessing  DSS.  Absence  of  these  routines 
or  the  presence  of  the  dummy  DSMSTR  does  not  have  any  computational  effect  on 
program  execution. 

3.0  ERROR  MESSAGES 

The  HEC-1  Users  Manual  (Section  11,  Table  11.1)  lists  error  messages  Which 
will  print,  along  with  an  explanation  of  the  message.  Some  errors  will  not 
cause  the  program  to  stop  execution,  so  the  user  should  always  check  the  output 
for  possible  errors  or  warnings.  The  array  dimensions  listed  in  Table  IV  are 
those  used  by  HEC  on  a  CDC  Cyber  865  computer,  and  correspond  to  the  official 
HEC-distributed  version  of  the  program.  Other  dimensions  may  be  used  as  noted 
in  the  following  section  on  Memory  Management.  If  the  dimensions  are  changed, 
the  dimension-error  checks  should  also  be  changed  in  the  source  code. 

The  computer  operating  system  may  also  print  error  messages.  When  an  error 
occurs,  the  user  should  first  ascertain  if  it  is  generated  by  HEC-1  or  by  the 
system.  If  it  is  generated  by  HEC-1,  i.e.,  in  the  format  given  in  the  Users 
Manual,  that  table  should  be  referred  to  and  the  indicated  actions  taken.  If 
the  error  is  system  generated,  the  computer  center  user  service  and/or  the  in- 
house  computer  systems  personnel  should  be  contacted  to  ascertain  the  meaning 
of  the  error.  These  errors  may  be  due  to  incorrectly  input  or  read  data, 
errors  in  HEC-1,  or  the  computer  system.  If  these  system  errors  cannot  be 
resolved  in-house  or  if  there  appears  to  be  an  error  in  the  HEC-1  program,  the 
HEC  should  be  contacted. 

4.0  COMPUTER  MEMORY  MANAGEMENT 

4.1  ChanRinx  Dimensions 

Computer  memory  requirements  can  be  reduced  by  decreasing  the  dimensions 
of  some  arrays.  Table  IV  identifies  the  arrays  whose  dimensions  have  the  most 
impact  on  core  storage.  Because  of  the  error  checks  on  dimensions  (see 
previous  section)  and  for  ease  of  making  changes,  the  dimensions  are  specified 
as  variables . 

To  change  a  dimension,  both  the  dimension-size  variable  and  the  arrays 
themselves  must  be  changed.  Table  IV  shows  both  the  dimension-size  variables 
and  the  arrays  in  which  the  particular  dimension  occurs.  The  labeled  common 
block  or  subroutines  in  which  the  array  occurs  is  also  identified.  Note  that 
many  of  the  dimensions  are  a  function  of  other  dimension  statements.  The 
amount  of  reduction  in  a  dimension  will  be  dependent  upon  the  options  and  the 
sizes  desired  by  the  user. 

Pay  special  attention  to  arrays  which  occur  in  EQUIVALENCE  statements. 

4 . 2  Deleting  Subroutines 

Computer  memory  requirements  can  also  be  reduced  by  eliminating  some 
subroutines.  The  subroutines  required  for  various  HEC-1  simulation  options 
are  shown  in  Table  V.  The  memory  reduction  obtained  from  deleting  subroutines 
may  not  be  significant  if  an  overlay  is  used. 
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Table  IV 


Array  Dimensions 


Variable 

Array 

Name/Definition  Size  Limit 

Arrays 

Location 

KQ 

Hydrograph  ordinates 

300 

DEWPT , EXCSR , EXCSS , ITLS 
KDAY , KHOUR , KMON , PRCP , 
PRCPR , Q , QO , RAINA , SNMT , 
SOL , STG , STR , TMPR , WIND 

COMMON/ORDTS/ 

KR 

Recording  precip  gages 

15 

PRCPR 

COMMON/ORDTS/ 

I  STAR 

COMMON/PRECIP/ 

KQH1 

PLAN/RATIO  hydrograph 
ordinates 

4800 

QH.QK 

COMMON/MULTQ/ 

KUHGQ 

Unit  graph  ordinates 

150 

QUNGR , QCLK 

COMMON/UNITGR/ 

KN 

Non-recording  gages 

70 

ANAPN , I STAN 

COMMON/PRECP/ 

KHN 

Non-recording  gage 
weighting 

10 

ISTN.WTN 

COMMON/PRECP/ 

KHR 

Recording  gage  weighting 

5 

ISTR.WTR 

COMMON/PRECP/ 

KZONE 

Snowmelt  zones 

10 

ANAP , ANDAY , AREA , CUML , 
SNO , SNOW 

COMMON /SNOWZ/ 

KSTM 

Depth-area  precip  array 

9 

AX , BX , CX , DX , STRM , TRDA 

COMMON/MULTPS/ 

KRTIO 

Flood/precip  RATIOS 

9 

RTIO 

COMMON /MULTPS/ 

FREQ , PFRQ 

COMMON/ECON2/ 

QPREP 

SUBROUTINE  ECONO 

LCM2 

Internal  Random  Access  24025 
Array 

SAVE 

COMMON/RAND/ 

XKQH  should  be  the  lesser  of  2*KRTI0»KQ  or  4800 
2LCM  =  5  •  ( KQH+5  ) 
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4 .3  Overlaying  Subroutines 


Computer  memory  requirements  can  also  be  reduced  by  overlaying  subroutines. 
A  suggested  overlay  structure  is  given  in  Table  VI.  That  overlay  structure 
was  successfully  used  to  implement  the  program  on  a  Harris  100  computer  with 
64k  (decimal)  words  of  central  memory. 

The  overlay  structure  details  which  subroutines  must  be  in  memory  at  the 
same  time  during  execution  by  a  root  path.  The  path  is  described  by  the  letter 
code  given  in  Table  VI.  For  example,  the  MAIN  ROOT  connects  with  high  level 
sub  roots  A,  B,  C,  D,  E,  F,  G,  or  H.  The  sub- roots  connect  with  lower  level 
roots.  For  example,  sub- root  B  connects  with  Bl,  B2 ,  B3 ,  B4,or  B5 .  This  root 
heirarchy  may  be  continued  to  lower  sub  root  levels.  For  example,  a  root  path 
may  consist  of  F,  F3 ,  F3B  and  either  F3B1 ,  F3B2 ,  F3B3,  F3B4 ,  or  F3B5 .  All  the 
subroutines  that  are  part  of  any  root  path  must  be  in  memory  at  the  same  time. 

5.0  PROGRAM  STRUCTURE 

The  HEC  1  computer  program  consists  of  a  main  routine  and  67  subroutines. 
The  interlinkage  of  the  subroutines  is  shown  in  Figure  1.  A  brief  description 
of  each  subroutine  is  given  in  Table  VII. 
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Table  V 


Required  Subroutines  for  HEC-1  Simulation  Options 


Simulation  options 

1  All  jobs 

2  Runoff  parameter  optimization 

3  Routing  parameter  optimization 

4  Stream  network 

5  Multiplan-multiratio 

6  Economics 

7  System  optimization 

8  Stream  network  with  dam  break 

9  Stream  nework  with  kinematic  wave 

10  Stream  network  with  DSS 

11  Economic  model  with  DSS 

X  required  subroutines 

*  subroutine  occasionally  used,  depends  on  input  option 
f  subroutines  needed  for  DSS  (Corps  users  only) 


Options 


1 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

Subroutine! 

HECl 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

ADDTIM 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

ADJUST 

1 

* 

* 

* 

* 

* 

* 

* 

* 

AKIMA 

1 

X 

X 

X 

AKIMAI 

1 

X 

X 

X 

BALNC 

1 

* 

BANNER 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

BASIN 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

BIT 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

CHNPRC 

1 

* 

* 

* 

* 

ft 

* 

* 

* 

CKPRNM+ 

1 

X 

X 

COPYIN 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

COST 

1 

X 

DAMAGE 

1 

X 

X 

X 

DAKBRK 

1 

X 

DATES 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

DIVERT 

1 

X 

X 

X 

X 

X 

X 

X 

X 

DSMSTR 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

EADSUM 

1 

X 

X 

X 

EAV 

1 

X 

X 

X 

ECONO 

1 

X 

X 

X 

EQUATN 

1 

* 

* 

* 

X 

X 

* 

* 

* 

EREXIT 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

FDKRUT 

1 

X 

FDKSRM 

1 

X 

FLOGRD 

1 

X 

FRMMTC 

1 

X 

GETDTA 

1 

X 

GRAPH 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

HYDLMT 

1 

* 

* 

* 

* 

X 

* 

* 

* 
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Table  V  (Continued) 


Simulation  options 

1  All  jobs 

2  Runoff  parameter  optimization 

3  Routing  parameter  optimization 

4  Stream  network 

5  Multiplan-multiratio 

6  Economics 

7  System  optimization 

8  Stream  network  with  dam  break 

9  Stream  nework  with  kinematic  wave 

10  Stream  network  with  DSS 

11  Economic  model  with  DSS 

X  required  subroutines 

*  subroutine  occasionally  used,  depends  on  input  option 
t  subroutines  needed  for  DSS  (Corps  users  only) 

Options 


1  1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

Subroutine  1 

IDTC 

1  x 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

INPUT 

1  X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

INTPL 

|  X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

INTPLQ 

1 

* 

* 

* 

* 

* 

* 

* 

A 

INVAR 

1 

X 

X 

KINOFF 

1 

X 

LOCALQ 

1 

* 

* 

* 

LOCAT 

1 

* 

* 

* 

* 

* 

* 

* 

A 

LOGLOG 

1 

* 

* 

* 

* 

* 

* 

* 

A 

LPDMG 

1 

X 

LTR 

1  X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

MANE  2 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

MISDTA 

1 

* 

* 

X 

X 

OPSTAT 

1 

X 

OPTIM 

1 

X 

X 

X 

OUTPUT 

1  x 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

OUTSUM 

1 

X 

X 

X 

X 

X 

X 

X 

X 

PLNSET 

1 

X 

X 

X 

X 

X 

PMPSAV 

1 

* 

* 

* 

X 

* 

* 

* 

A 

PREVU 

1 

X 

X 

X 

X 

X 

X 

X 

X 

PRNTQ 

1  x 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

PRT 

1 

X 

X 

X 

PUMP 

1 

* 

* 

* 

X 

* 

* 

A 

A 

PSUMRY 

1 

* 

* 

* 

X 

A 

* 

A 

A 

QG3 

1 

X 

X 

X 

RDHDPR+ 

1 

X 

X 

RDTIMS+ 

1 

X 

X 

RDZWEC+ 

1 

X 

RDZWTS+ 

1 

X 

X 

READBZ+ 

1 

X 

X 

READQ 

1 

* 

A 

* 

* 

* 

* 

A 

A 

READZR+ 

1 

X 

X 

RECESS 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 
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Simulation  options 


Table  V  (Continued) 


1  All  jobs 

2  Runoff  parameter  optimization 

3  Routing  parameter  optimization 

4  Stream  network 

5  Multiplan-multiratio 

6  Economics 

7  System  optimization 

8  Stream  network  with  dam  break 

9  Stream  nework  with  kinematic  wave 

10  Stream  network  with  DSS 

11  Economic  model  with  DSS 

X  required  subroutines 

*  subroutine  occasionally  used,  depends  on  input  option 
+  subroutines  needed  for  DSS  (Corps  users  only) 


Options 


1  1 
l 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

Subroutine! 

REDARY 

1  A 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

RESOUT 

1 

X 

X 

X 

X 

X 

X 

X 

X 

RESQ 

1 

* 

* 

* 

X 

X 

* 

* 

* 

RESVRT 

1 

* 

* 

* 

X 

X 

* 

* 

* 

ROFGRD 

1 

X 

ROUTE 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

RUNOF 

1 

X 

X 

X 

X 

X 

X 

X 

X 

X 

SAVDTA 

1 

X 

SETBIT 

1  A 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

SETOPT 

1 

X 

SOLVEQ 

i 

* 

* 

* 

X 

X 

* 

* 

* 

SPIRAT 

! 

* 

* 

* 

* 

* 

* 

* 

* 

STACK 

1 

X 

X 

X 

X 

X 

X 

X 

X 

STAGE 

1 

* 

* 

* 

* 

* 

* 

STFN 

1 

* 

* 

* 

A 

X 

* 

* 

* 

STORAG 

1 

* 

* 

* 

X 

X 

* 

* 

* 

STOUT 

1 

* 

* 

* 

* 

X 

* 

* 

* 

STPRTC+ 

1 

X 

X 

SUMOP 

1 

X 

SUMRY 

1  A 

X 

X 

X 

X 

X 

A 

X 

X 

A 

X 

SYSOPT 

1 

X 

TCOMPR 

1  A 

A 

X 

A 

X 

X 

X 

X 

X 

A 

X 

TIMAVG 

1  A 

A 

A 

A 

X 

X 

X 

X 

X 

A 

X 

TIMEIN 

j  A 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

TOMTRC 

1 

X 

TP40 

1 

* 

* 

* 

* 

* 

* 

* 

* 

TSUB 

1  A 

X 

A 

A 

X 

X 

X 

X 

X 

A 

X 

TX2IN 

|  A 

X 

X 

A 

X 

X 

X 

X 

X 

X 

X 

USRDO 

1 

* 

* 

* 

* 

* 

* 

WRITQ 

1 

* 

* 

* 

* 

* 

* 

WRTQDS+ 

1 

X 

X 

WRTQFR+ 

1 

X 

ZERO 

1  A 

A 

X 

A 

X 

X 

X 

X 

X 

X 

X 

ZRTYPE+ 

1 

- 

X 

X 
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Table  VI 

Overlay  Structure 


Main  Root 
HEC1 
DSMSTR 
EREXIT 
OPTIM 
SAVDTA 
SETOPT 
SYSOPT 

USRDO 

| (X) BANNER  COPY IN 

PREVU 

PRT 

ZERO  | 

i  d) 

|  BIT 
j  IDTC 
j  INPUT 
j  LTR 
j  REDARY 
j  SERBIT 
|  TCOMPR 
|  TSUB 

1 

|  (Bl)  DATES 

|  0)  TIMEIN 

ADDTIM 

TX2IN 

1 

I 

|  0)  ADJUST 

BASIN 

LOCAT 

LOGLOG 

TP40 

|  0  INTPLQ 

PRNTQ 

READQ 

SUMRY 

TIMAVG 

|  (£5)  MISDTA 
j  W  READBZ 

CKPRNM 

READZR 

RDHDPR 

STPRTC 

RDTIMS 

ZRTYPE 

RDZWTS 

(c)bit  OUTPUT  I 


BIT  RESOUT  DAMBRK  |  D1  HYDLMT  STOUT  | 

(5)  COST  RESQ  GRAPH  j  ===========_===^  _-== 

INTPL  STORAG  j  D2  SPIRAT  | 

©BALNC  PRNTQ  STACK  SUMRY  TIMAVG  | 


821=33= 


PMPSAV 

TIMAVG 


(6)AKIMA  AKIMAI  BIT  COST  DAMAGE  EADSUM  EAV  | 

U ECONO  INTPL  LPDMG  QG3  SETBIT  RDZWEC  WRTQFR  | 

Sa33S8=SCSS8SSS388S3a8SSaS9«S3S3SS3S38SffSS9SBa3SS3aSS3H8£S3a 

^h)getdta  I 


PSUMRY  | 
OUTSUM  j 


Flgur  •  1  INTERCONNCCTION 


Table  VII 


HEC1 

ADDTIM 

ADJUST 

AKIMA 

AKIMAI 

BALNC 

BANNER 

BASIN 

BIT 

CHNPRC 

CKPRNM 

COPYIN 

COST 

DAMAGE 

DAMBRK 

DATES 

DIVERT 

DSMSTR 

EADSUM 

EAV 

ECONO 

EQUATN 

EREXIT 


Description  of  HEC-1  Subroutines 

initialization,  control  program  flow 

adds  one-time  interval  to  current  time 

computes  depth-area  adjustment  for  hypothetical  storm 

finds  coefficients  for  cubic  spline 

cubic  spline  interpolation 

computes  balanced  hydrograph 

prints  program  banner 

computes  basin-average  rainfall;  standard  project  storm;  provable 
maximum  storm;  weighting  of  gaged  rainfall 

tests  a  bit  in  an  integer  variable 

computes  channel  loss 

checks  parameter  name  for  DSS  variables  against  known  parameters 

copies  input  data  to  working  file;  adds  identification  of  next  record 
to  beginning  of  each  record;  converts  data  from  free  format  to  fixed 
format 

computes  cost  of  a  project  element 
computes  damage-frequency  relation 

computes  dam-breach  size;  sets  time  step  for  breach  calculation; 
prints  dam-break  summary 

computes  date  and  time  for  each  hydrograph  ordinate 
diverts  flows;  retrieves  diversion  hydrographs 
controls  use  of  DSS 

prints  summary  of  expected  annual  damages 
computes  expected  annual  damage 
reads  data  for  expected  annual  dai.tage  calculation 
equations  which  require  iterative  solution 
prints  number  of  errors;  exits  program  on  error 


Table  VII  (Continued) 


FDKRUT 

FDKSRM 

FLOGRD 

FRMHTC 

GETDTA 

GRAPH 

HYDLMT 

IDTC 

INPUT 

INTPL 

INTPLQ 

INVAR 

KINOFF 

LOCALQ 

LOCAT 

LOGLOG 

LPDMG 

LTR 

MANE  2 

MISDTA 

OPSTAT 

OPTIM 

OUTPUT 

OUTSUM 

PLNSET 


finite  difference  solution  of  kinematic  wave  routing 
routes  flows  using  kinematic  wave  method 

computes  distance  step  for  kinematic  wave  routing  of  runoff  hydrograph 
in  collector  and  main  channels 

converts  excess  from  metric  to  English  units 

gets  data  to  compute  plan  2  for  flood-contol  system  optimization 
makes  printer  plots  of  hydrographs 

computes  hydraulic  characteristics  from  cross-section  data 

compares  next  record  identification  code  with  current  record 
identification  code 

reads  input  data  from  working  file 

linear  interpolation 

interpolates  a  hydrograph  from  one  interval  to  another 

initializes  variables  for  runoff  parameter  or  routing  parameter 
optimization 

transforms  rainfall  excess  to  t/unofc  using  kinematic  wave 
computes  local  flow  hydrograph 

generates  position  in  triangular  distribution  for  hypothetical  storm 
log-log  interpolation 

computes  flow/ stage-damage  relation  for  current  local  protection  level 

finds  position  of  record  identification  code  in  array  of  valid  letters 

directs  program  to  desired  operation;  combines  hydrographs 

tests  a  variable  for  a  value  indicating  missing  data 

computes  statistics  for  comparing  two  hydrographs 

controls  variation  of  parameters  being  optimized 

prints  input  data 

prints  final  hydrograph  summary 

sets  variable  for  computing  remaining  plans  after  all  input  data  has 
been  read  for  a  station 
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Table  VII  (Continued) 


PMPSAV 
PREVU 

PRNTQ 
PRT 

PSUMRY  - 

PUMP 

QG3 

RDHDPR 

RDTIMS  - 

RDZWEC 

RDZWTS 

READBZ 

READQ 

READZR 

RECESN 

REDARY 

RESOUT 

RESQ 

RESVRT 

ROFGRD 

ROUTE 

RUNOF 

SAVDTA 

SETBIT 


saves  pumpflow  hydrographs 

prints  schematic  diagram  of  stream  system;  sets  up  table  for  user- 
defined  output 

prints  hydrographs 

copies  summary  table  produced  by  subroutines  DAMBRK  and  SUMOP  from 
scratch  files  to  printer 

saves  pumpflow  values  for  summary  printout 

computes  pump  flow 

three  point  gaussian  quadrature 

reads  HYDPAR  data  from  DSS 

reads  time  series  data  from  DSS 

-  reads  ZW  record  for  f low-frequency  data 

-  reads  ZW  record  for  hydrologic  time  series 

-  reads  BZ  record 

-  reads  a  hydrograph  from  a  file 

-  reads  ZR  record 

-  adds  base  flow  to  direct  runoff  hydrograph 

-  reads  table  data  (elevation,  flow,  storage,  etc.) 

-  computes  spillway  flow  from  weir  coefficients;  computes  storage- 
outflow  table  from  storage-elevation  and  outflow-elevation  tables 

-  computes  flow  from  reservoir  for  given  elevation 

-  routes  flow  through  reservoir  using  level-pool  routing 

-  computes  time  and  distance  step  size  for  kinematic  wave  routing  of 
rainfall  excess 

-  routes  flows  using  Muskingum,  modified  Puls,  or  straddle-stagger 
methods;  optimizes  Muskingum  K  and  X 

-  computes  unit  hydrograph;  precipitation  excess;  transforms  excess  to 
direct  runoff 

-  saves  data  for  computing  plan  2  for  flood-control  system  optimization 

-  sets  a  bit  in  an  integer  variable 
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Table  VII  (Continued) 


SETOPT 

SOLVEQ 

SPIRAT 

STACK 

STAGE 

STFN 

STORAG 

STOUT 

STPRTC 

SUMOP 

SUMRY 

SYSOPT 

TCOMPR 

TIMAVG 

TIMEIN 

TOMTRC 

TP40 

TSUB 

TX2IN 

USRDO 

WRITQ 

WRTQDS 

WRTQFR 

ZERO 

ZRTYPE 


sets  variable  to  be  optimized 

iterative  solution  algorithm  using  combination  of  secant  method  and 
interval  halving 

computes  flow  from  ogee  spillway 
stores  and  retrieves  hydrographs 
converts  flows  to  stages 

computes  change  in  storage  plus  inflow  minus  outflow  (continuity)  for 
given  elevation  for  reservoir  routing 

computes  reservoir  storage  from  surface  area  and  elevation  data 

computes  storage-outflow-elevation  table  for  modified  Puls  routing 
from  eight-point  cross-section 

sets  part  C  of  DSS  pathname;  saves  parameter  name 
prints  summary  of  runoff  parameter  optimization 
computes  summary  of  a  hydrograph  computation 
controls  optimization  of  flood-control  system  components 
compares  two  times  in  internal  form 
computes  average  values  for  a  time  series 

reads  time  series  data;  interpolates  to  computation  interval 
converts  runoff  hydrograph  from  English  to  metric  units 
computes  rainfall  distribution  for  hypothetical  storm 
subtracts  two  times  in  internal  form 

converts  external  time  (day,  month,  year,  hour)  to  internal  time 
(year,  day,  minute) 

prints  user-defined  output  tables 

writes  a  hydrograph  to  a  file 

writes  a  hydrograph  to  DSS 

writes  flow-frequency  data  to  DSS 

initializes  variables  for  a  new  job 

gets  replacement  record  type  from  a  ZR  record 
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